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ABSTRACT 

Z sources are bright neutron-star X-ray binaries, accreting at around the Eddington limit. We 
analyze the 68 RXTE observations (^270 ks) of Sco-like Z source GX 17+2 made between 1999 
October 3-12, covering a complete Z track. We create and fit color-resolved spectra with a model 
consisting of a thermal multicolor disk, a single-temperature-blackbody boundary layer and a weak 
Comptonized component. We find that, similar to what was observed for XTE J 170 1-462 in its Sco- 
like Z phase, the branches of GX 17+2 can be explained by three processes operating at a constant 
accretion rate M into the disk: increase of Comptonization up the horizontal branch, transition from 
a standard thin disk to a slim disk up the normal branch, and temporary fast decrease of the inner 
disk radius up the flaring branch. We also model the Comptonization in an empirically self-consistent 
way, with its seed photons tied to the thermal disk component and corrected for to recover the pre- 
Comptonized thermal disk emission. This allows us to show a constant M along the entire Z track 
based on the thermal disk component. We also measure the upper kHz QPO frequency and find 
it to depend on the apparent inner disk radius i?i n (prior to Compton scattering) approximately 

as frequency cx R in , supporting the identification of it as the Keplerian frequency at R m . The 
horizontal branch oscillation is probably related to the dynamics in the inner disk as well, as both its 
frequency and i?; n vary significantly on the horizontal branch but become relatively constant on the 
normal branch. 

Subject headings: accretion, accretion disks — starts: individual (GX 17+2) — stars: neutron - 
X-rays: binaries — X-rays: bursts — X-ray: stars 



1. INTRODUCTION 

^ ' Based on the timing and spectral properties, six of 
. the persistently bright neutron-star (NS) low-mass X- 
• ray binaries (LMXBs) were classified more than two 
*~ 1 ' decades ago as Z sources, named after the spectral evo- 
lution patterns they trace out in X-ray color-color di- 
r~^. ■ agrams (CDs) or hardn ess-i ntensity diagrams (HIDs; 
Q 1 IHasinger fc van der Klii 19891: Ivan der Klisll2006l) . These 
(NJ ! sources are Sco X-l, GX 17+2, GX 349+2, GX 340+0, 
| ■ GX 5-1, and Cyg X-2. The upper, diagonal and lower 
branches of their Z-shaped tracks are called horizontal, 
. £^ | normal, and flaring branches (HB/NB/FB), respectively. 
K> . Based on the shape and orientation of the Z tracks, these 
Z sources were further divided into two subgroups, with 
d ' the first three called "S co-like" Z sources and the latter 
three called "Cyg-like" (IKuulkers et al.lll99l . The Sco- 
like Z sources have a more vertically oriented HB and a 
stronger FB than the Cyg-like types. The Z tracks them- 
selves can also move and change shapes in the CDs/HIDs 
(secular changes), most substantially in the case of Cyg 
X-2. 

Recently, XTE J 170 1-462 became the first X-ray tran- 
sient identified as a Z-source, and studies of this source 
have significantly improved our understanding of the 
spectral evolution in Z sources and their relation to atoll 
sources, another subclass of NS LMXBs with lower X-ray 
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luminosity (Lx) than Z sources. XTE J1701-462 experi- 
enced a long outburst in 2006-2007, and it showed suc- 
cessive characteristics of Cyg-like Z, Sco-li ke Z and atoll 
sources dur i ng the decay of the outburst (jHoman et al.l 
120071 l2ljll ILin et al.l [20091 LRH09 hereafter). During 
these secular changes in the Z tracks, the upper (HB/NB) 
and lower (NB/FB) vertices each evolved along a spe- 
cific line in the HID (LRH09). All three Z branches were 
present when the source was bright. During the decay, 
the HB, NB, and FB successively disappeared, and the 
lower vertex finally transitioned into the atoll soft state. 

Using an X-ray spectral mo del that was su ccessfully 
applied to two atoll transients (|Lin et al.ll2007fl . LRH09 
showed that it is most likely the mass accretion rate (M) 
into the disk that drives the secular changes of Z tracks 
and the transitions from Z to atoll types. While the inner 
disk radius remained constant in the soft state of the 
atoll stage from the same spectral model, the inner disk 
radius in the Z stage increased with luminosity, which 
was interpreted as an effect of the local Eddington limit. 
On the other hand, the motion along a Z branch on short 
timescales appeared to operate at roughly constant M, 
at least for Sco-like Z tracks (LRH09). Furthermore, the 
three Z branches were linked to different mechanisms by 
LRH09. The source ascends the HB as Comptonization 
of the disk emission increases. The apparent luminosity 
of the boundary layer increases along the NB from the 
lower to the upper vertices, which can be explained as 
a transition from a geometrically thin disk to a thick 
disk that is expected to admit an advective component 
of mass flow to the NS. The FB is traced out as the inner 
disk radius temporarily decreases toward the value seen 
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Fig. 1. — RXTE ASM one-day-averaged light curve of GX 17+2 spanning ~14 years. The narrow gray region marks the interval UTC 
1999-10-03.1-12.3, during which the PCA observations analyzed in the paper were obtained. 



in the atoll soft state, presumably the innermost stable 
circular orbit (ISCO). 

In this work we study GX 17+2. Figurc[T]plots its long- 
term light cu rve from the RXTE All-Sky Monitor (ASM; 
ILevine et all l996). The source shows a roughly constant 
level (~45 counts/s), with frequent flares on top, suggest- 
ing little change of the overall properties of this source 
over a long ti mescale. In fact, it shows very small sec - 
ular changes (jWiinands et al.1119971; IHoman et al.ll2002[) . 
in contrast to XTE J1701-462. 

X-ray Spectral s tudies of GX 17+2 have been re- 
ported previously dDi Salvo et al.l l2000t iFarinelli et al.l 



2(11).-): IMiriiari et all 120071) . IDi Salvo et al. I (12000D fit- 
ted the BeppoSAX spectra from the HB and NB with a 
single-temperature blackbody (BB) plus a Comptonized 
component for the continuum spectra. They found a 
hard tail on the HB, which was fit by a power law 
(PL) and contributed -8% of the 0.1-200 keV source 
flux. This component gradually faded as the source 
moved toward th e NB, where it wa s no longer de- 
tectab le (see also Farinclli et al.l [20051 ). iMigliari et al.l 
(|2007t ) carried out simultaneous radio and X-ray obser- 
vations of GX 17+2 and found that a positive correla- 
tion between the radio flux density and the flux of the 
hard tail . Detailed tim i ng an a lyses of GX 17+2 wer e 
made bvlKuulkers et all (|1997l ); IWiTnands et al.l (|1997t ); 
IHoman et al.l (|2002l ). with the evolution of various types 
of quasi-periodic oscillations (QPOs) along the Z track 
obtained. 

Here we carry out the spectral modeling of GX 17+2, 
in order to determine whether its spectral evolution is 
similar to that of the Sco-like Z stage in XTE J1701- 
462 (LRH09). As its secular changes are small, we 
just concentrate on a 9-day period of intensive RXTE 
pointed observations in 1999. Different from previous 
X-ray spectral modeling of this source, we use a low- 
Comptonization model dominated by a thermal disk and 
a thermal boundary layer component, because of its suc- 
cesses with atoll sources and XTE J1701-462 (|Lin et al.l 
120071 120101 LRH09). One main advan tage of this m odel 
over other commonly used models (see lBarretll200lL for a 
review) is that it infers the disk in the soft state of atoll 
sources to behave approximately as L oc T 4 , which is of- 
ten seen in the thermal state of black hole X-ray binaries 



as well (|Lin et al.l 12007). We also compare our spectral 
fit results with the evolution of kHz QPOs and HB os- 
cillations (HBOs) in order to obtain hints on the origins 
of these QPOs. GX 17+2 is more suitable than XTE 
J1701-462 for such a study, since such Q POs were not 
detected as frequently in the latter source ([Homan et al.l 
l2llltlSannalt^l2mol:rBarret et al.ll2mT[ ). In Section^ 
we describe the reduction of the data and the procedure 
by which we create our spectra. The CDs and HIDs are 
also presented in this section. We describe the spectral 
models in Section [3] The spectral fit results and cor- 
relations with fast varibility are given in Section 2J In 
Section [51 we discuss the mass accretion rate, the physi- 
cal interpretation of spectral evolution, and the possible 
origins of kHz QPOs and HBOs. Finally we present our 
conclusions. 

2. DATA ANALYSES AND COLOR-COLOR DIAGRAMS 

We analyzed 68 RXTE observations of GX 17+2 made 
between UTC 1999 October 3.1-12.3 (i.e., MJD 51454.1- 
51463.3), using the same standard criteria to filter the 
data (e.g. , removal of fiv e type I X-ray bursts) as de- 
scribed in lLin et al.l ([20071 ). FTOOLS 6.9 was used. The 
source can move on timescales as short as minutes in 
the CD/HID, especially on the FB. To track the evolu- 
tion along the Z track but also to gain enough statis- 
tics for spectral modeling and timing analysis, data are 
normally split into short exposures and then rebinned 
based on the i r pos itions in the CD/HID (e.g., LRH09; 
IHoman et al.l [2002) . We also followed this procedure 
here. We calculated X-ray colors using the Proportional 
Counter U nit (PCU) 2 of the Proportional Counter Ar- 
ray fPCA : Uahoda et all 119961) . in the same way as de- 
scribed in lLin et al.l (2002& . Th e soft color (SC and the 
hard color (HC) are defined as the ratios of the Crab- 
Nebula-normalized count rates in the (3.6-5.0)/(2.2-3.6) 
keV bands and the (8.6-18.0)/(5.0-8.6) keV bands, re- 
spectively. We used spectra with exposures of 32 s from 
the "standard 2" data to calculate the colors and con- 
struct the CD and HID, which are shown in the upper 
panels of Figure [5] The HB, NB, and FB are denoted by 
the blue square, green triangle, and red cross symbols, 
respectively, and this convention is used for all other fig- 
ures in this paper. 
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Fig. 2. — The color-color and hardness-intensity diagrams. The HB, NB and FB are marked by blue squares, green triangles, red crosses, 
respectively. Upper panels: 32 s spectra are used, and the splines that are used for the Sz parametrization are shown in the HID. Lower 
panels: Sz-resolved spectra are used, and the splines are repeated here for reference. 
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Fig. 3. — RXTE PCA 32-s light curves in two energy bands. The typical error bars are smaller than the symbol size. The meanings of 
the symbols are the same as Figure [2] 



Only a single Z track is seen in these diagrams. The 
lower parts of the NB and FB have substantial overlap. 
From the 2.2-3.6 keV and 8.6-18.0 keV light curves in 
Figure 02 we see that the source experienced frequent 
flares between MJD 51461.2 and 51463.3. We identified 
the data within this interval as the FB, as supported 
by their positions in the CD and HID in Figure [2j The 
transition from the HB to the NB is quite smooth. Their 
boundary is defined to be HC ~ 0.45 (strictly Sz = 1, see 
below). We note that the branches are hard to separate 



completely and that the above identification is approxi- 
mate. 

The data points are often assigned a r ank number Sz 
to track their positions along the Z track (fl asinger et al.l 
119901: IHertz et alJU992T ). This can be done by creating 
splines along the Z track. Considering that the scatter- 
ing is large in the CD (mostly due to large statistical 
uncertainties of SCs), we created the splines in the HID, 
which are shown in the top right panel of Figure [2] Due 
to large overlap of lower parts of the FB and the NB, we 
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created a spline for the FB (the FB spline) and another 
one for the HB and NB (the HB+NB spline). We picked 
the normal points of the splines by hand, with one of 
them shared by the two splines. The Sz value of this 
common normal point is set to be 2.0 (the lower vertex), 
while the point with HC=0.45 in the HB+NB spline has 
Sz = 1.0 (defined as the upper vertex and the bound- 
ary the HB and NB). Sz values at other positions in the 
splines are determined based on their distances to these 
two points along the splines. Considering the different 
units of HC and intensity, we divide them by a charac- 
teristic number before calculating the distances, 0.45 for 
the HC and 2000 cts/s/PCU for the intensity. Finally, 
the values of Sz for data points on the FB and on the 
HB and NB are obtained by projecting them onto the 
F B and the HB+NB s plines, respectively. We note that 
in IHoman et al.l (|2002f) a single spline was used and was 
applied in the CD. This is due to wider energy bands 
that they used to define the SCs. We chose t o kee p our 
definition of colors as adopted in iLin et ail (|2007|) and 
LRH09 for easy comparison with these studies. 

We then combined 32 s spectra from PCU 2 based on 
selections in Sz to create spectra with longer exposures 
for our spectral fits. The ranges of Sz and the exposures 
for these Sz-resolved spectra are given in Tables [TH2J 
Their locations in the CD and HID are shown in the 
lower panels of Figure [2] We also created spectra of 
PCU and those o f the High Energy X-ra y Timing Ex- 
periment (HEXTE; [Rothschild et al.ll998D . which match 
the spectra of PCU 2 in time. PCUs and 2 were the 
only units that were operating during all of our observa- 
tions. For the HEXTE, we used both Clusters A and B. 
We applied dead time corrections to both the PC A and 
HEXTE data. 

We also created power density spectra (PDS), which 
match the Sz-resolved spectra in time, to search for the 
kHz QPOs and HBOs and study their relations with 
our spectral fit paramet ers. As we us e d a di fferent way 
of data reduction from IHoman et al.l |2002), our PDS 
do not match theirs in time. We used the PCA event 
files to create 5.8-40 keV light curves (from all avail- 
able PCUs) in a time bin of 2 -12 s. We then calculated 
Sz-resolved PDS by averaging the individual PDS com- 
puted for the appropriate set of 32-s time intervals. We 
excluded events below 5.8 keV to increase the d etection 
significance of QPOs, as in IHoman et al.l (|2002l) . We re- 
fer to their paper for more information on the data modes 
of the even t files. We then fitte d the PDS in a very sim- 
ilar way as IHoman et al.l (|2002l ) to search for kHz QPOs 
and HBOs. Basically the 200-2048 Hz PDS were fit- 
ted with one or two Lorentzians for the kHz QPOs plus 
a dead-time-modified Poisson level, and the 0.03125-256 
Hz PDS were fitted with a PL, a cutoff PL (CPL) and/or 
a Lorentzian for various broad noise components and 
Lorentzians for QPOs, including HBOs, plus a constant 
for the Poisson level. We also calculated the rms within 
the 0.1-10 Hz frequency band from these PDS, in order 
to study its correlation with the Comptonization fraction 
inferred from the spectral parameters. 

As the PCA does not cover the low energy band be- 
low 2.5 keV and has only a modest energy resolution, we 
also analyzed two Suzaku observations (sequence num- 
bers 402050010 and 402050020), using the same data 
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Fig. 4. — The fit residuls in terms of sigma for the spectrum 
with Sz =0.32-0.40, using various models for the Comptonization 
(see text for details). This spectrum is used considering its large 
Comptonization fraction and exposure. (A): a PL, x2=l-13; (B) 
a CPL, i.e., Model 1, x 2=0.80; (C) an nthComp, i.e., Model 2, 
X?=0.86. The black, red, green, and blue data points are for the 
PCU2, PCU0, HEXTE Clusters A and B, respectively. 

reduction procedure as in ILin et al.l (|2010D . to help to 
estimate the absorption and the Fe line energy for our 
spectral fits to the RXTE data. The calibration files of 
2010 June and Su zaku FTOOLS versi on 15 were used. 
Spectra from XIS (|Kovama et al.ll2~007l ) and 3 were ex- 
tracted, with the central region of a radius of 55 pixels 
excluded to reduce the event pileup effect. XIS 1 was 
not used as it was in the full window mode, resulting 
in serious event pileup. We also used t he spectra of the 
PIN diodes of the Hard X-ray Detector ([Takahash i et al.1 
l200l . 

3. SPECTRAL MODELING 

The Sz-resolved PCA (PCU0 and PCU2) and HEXTE 
(Clusters A and B) spectra were fitted jointly, with their 
relative normalizations allowed to float. We used the 
energy range 2.9-40.0 keV for the PCA spectra and 20.0- 
60.0 keV for the HEXTE spectra. We assumed a model 
systematic error of 0.5%, as recommended by the PCA 
team. Both PCA and HEXTE data were binned to have 
at least 40 counts per bin. All error bars of spectral 
fit results quoted are at a 90%-confidence level, unless 
indicated otherwise. 

We fitted the spectra with a model consisting of a 
BB (bbodyrad in XSPEC), used to describe the bound- 
ary layer, a multicolor disk blackbody (MCD; diskbb in 
XSPEC), and a weak Comptonized component. It also 
includes a Fe K Gaussian line, a Fe K absorption edge, 
and an interstellar medium absorption, described by the 
gaussian, edge, and wabs models in XSPE C, respectively. 
The a bsorption edge was included in iDi Salvo et al.l 
(2000) in their fits to the BeppoSAX spectra, and we 
found that including it in our model also improved our 
fits significantly (% 2 decreases by about 37 on average 
for each Sz-resolved spectrum, with 130 degrees of free- 
dom typically). To reduce the scattering of our results, 
we used fixed values of the hydrogen column density 
TVh = 2.35xl0 22 cm~ 2 , the Gaussian line central en- 
ergy E ga = 6.64 keV, and the edge energy i? c d gc = 8.93 
keV (derived below and at the end of this section) for 
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(a) Model 1: MCD+BB+CPL 
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(b) Model 2: MCD+BB+nthComp 



Fig. 5. — The unfolded spectra and fit residuals for Models 1 and 2 at different positions in the Z track. For the unfolded spectra, 
only spectra of PCU2 and HEXTE Cluster A are shown for clarity. The total model fit is shown as a black solid line, the (unscattered) 
MCD component as a red dotted line, the BB component as a blue dashed line, the Comptonized component (CPL/nthComp) as a green 
dot-dashed line, and the Fe line as a cyan triple-dot-dashed line. For the fit residuals, the black, red, green, and blue data points are for 
the PCU2, PCU0, HEXTE Cluster A and B, respectively. 



all fits. We explored several descriptions of the Comp- 
tonized component. It is significant on the HB, and un- 
derstanding this component is important toward our un- 
derstanding of the HB, the HBOs, and the kHz QPOs. 

The use of three spectral components to fit X-ray spec- 
tra of acc reting NSs in the soft state must be done care- 
fully (see iLin et all I2007T ) . The Comptonized and the 
thermal components can become degenerate when the 
Comptonized component is steep (e.g., the photon index 
larger than 2.5 when the Comptonization is modeled as 
a PL). Here, we first used a CPL {cutoffpl in XSPEC) to 
describe the Comptonized component. The whole model 
is wabs(diskbb + bbodyrad + cutoffpl + gaussian)edge us- 
ing the XSPEC terminology, and this is our first model 
(Model 1, or model MCD+BB+CPL, hereafter). Using 
this model, the Comptonized component turned out to 
be weak compared with the thermal components except 
at the top of the HB. The photon index Tcpl and the e- 
folding energy Eqpl were generally not well constrained 
from the fits, and they showed no clear sign of variation 
along the Z track. To have a better constrained estimate 
of them, we fitted all the HB spectra simultaneously with 
these two parameters tied to be the same and obtained 
r CPL = 1.40+0.14 and E C pl = 9.9+1.0 keV. These val- 
ues were fixed in the final fits to all spectra. We then 
obtained £ e dge — 8.93 ± 0.03 keV from the simultaneous 
fit to the spectra from all branches. Although fits with 
the Comptonized component described by a PL can also 
have the xt values around one, such fits have systematic 
residuals above 20 keV on the HB (Figure HJ, which are 
significantly reduced using a CPL. The x 2 decrease is 
51.7, from 173.5 (d.o.f=153) to 121.8, for the S z =0.32- 
0.40 spec trum. Applying the posterior predictive p- va lue 
method (|Hurkett et al.1 l2008t IProtassov et all 120021 ) to 



this spectrum, we found that the reduction of x 2 from 
the introduction of Ecpl is less than 15 for all 10 4 spec- 
tra that we simulated. Thus the exponential rollover is 
required for the Comptonized component at a confidence 
level above 99. 99% (at least for the 5 z =0.32-0.40 spec- 
trum). Figure 5(a) shows the unfolded spectra and fit 



residuals at four representative positions of the Z track 
(the top of the HB, the upper vertex, the lower vertex, 
and the top of the FB) using Model 1. 

Next, we attempted to describe the Comptonization 
in a self-consistent way. Considering that the CPL lumi- 
nosity variation along the HB is strongly anti-correlated 
with the disk luminosity from Model 1 (Section |4} , we 
assumed the scenario that there is a corona above the 
disk so that some photons from the thermal disk emis- 
sion are scattered by the hot electrons in it and turned 
into Comptonization emission. Such a picture has been 
suggested in some studies of black hole X-ray binaries 
and can be modeled with the SIMPL model (in XSPEC, 
ISteiner et all I2009bl lah . SIMPL is an empirical convo- 
lution model of Comptonization in which a fraction of 
the photons from an input seed spectrum are scattered 
into a power-law component with the rest unscattered 
and observed directly. While applying the SIMPL model 
with the MCD as the input seed spectrum to our spectra 
at the top of the HB where the Comptonization is the 
strongest, we found some systematic fit residuals above 
20 keV, with the data falling below the model prediction 
at high energies (similar to case A in Figure |4] using a 
PL to describe the Comptonization). A possible expla- 
nation for this is that there is a high-energy cut-off in 
the Comptonization emission (i.e., the corona tempera- 
ture is close to or within the energy range of our data), 
which was not included in the SIMPL model for simplic- 
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ity (|Steiner et a l. 2009b). To account for this, one option 
is to modify the SIMPL model to include a high-energy 
cut-off for the Comptonized component by multiplying 
it by, e.g., the highecut model in XSPEC. An alterna- 
tive option is to replace SIMPL with a Comptonization 
model with the corona temperature as a parameter, such 
as the nthComp model (in XSPEC. iZvcki et "all 119991: 
iZdziarski et"allll996t iLightman fc Zdziarskilll987n . The 
nthComp model approximates the Comptonization by 
solving the Kompaneets equation with a relativistic cor- 
rection to energy transfer between photons and electrons. 
The input seed photons can be a BB or a MCD in spectral 
shape, and we assumed the latter throughout the paper. 
The key parameters of this model include the asymptotic 
power- law photon index r n thComp, the corona electron 
temperature fcT 0jn thComp, the seed photon temperature, 
and the normalization A„ t hComp- 

In this paper, we present the results with the Comp- 
tonized component modeled with nthComp. We note 
that similar conclusions can also be drawn using the 
SIMPL model modified to include the high-energy cut- 
off. Then our second model (Model 2 or model 
MCD+BB+nthComp, hereafter) is: wabs(diskbb + 
bbodyrad+nthcomp + gaussian)edge, with the seed pho- 
ton temperature tied to the temperature of the MCD 
component /cImcd- With Model 2, our aim is to take 
into account the photons Compton scattered so as to 
track the behavior of the thermal disk emission prior 
to Compton scattering (denoted as MCDPS hereafter), 
as done in SIMPL. Therefore, we calculated the photon 
flux of the MCDPS by adding the photon flux of nth- 
Comp (the scattered part of the MCDPS) to that of the 
MCD component (the unscattered part) and then ob- 
tained the corresponding energy flux and normalization 
of the MCDPS by increasing those of the MCD com- 
ponent proportionally to the photon flux (the MCDPS 
temperature ^Tmcdps is the same as A;Tmcd)- 

To have a better constrained estimate of the pa- 
rameters of the Comptonized component, we fitted the 
HB spectra simultaneously with their P n thComp and 
kT e< nthComp tied to be the same, as we did for Model 1. 
We found that comparably good fits can be obtained 
using a large range of r n thCom P , i.e., 2.0-2.7, with 
kT e< nthComp of 5.6-7.3 keV correspondingly. This uncer- 
tainty is coupled to the uncertainty in the fraction of the 
thermal disk emission scattered into the Comptonization 
emission, with a larger scattering fraction inferred using 
a larger value of r n thComp- However, the parameters (the 
temperature and the normalization) of the MCDPS and 
the BB are much less affected. Allowing a high value 
of PnthComp might have the risk of large interference be- 
tween the thermal (especially the disk) and the Comp- 
tonized components, as simila rly obsery e d in t he fits to 
atoll-source X-ray spectra in iLin et all (|2007l ). In the 
end, we chose to present results using TnthComp = 2.3. 
Then the fit to all the HB spectra simultaneously gave 
kT c 

, nthComp 

= 6.04 ±0.23. The above values of r n thCom P 
and fcT ejn thComp were fixed in the final fits to all spe ctra. 



Sample fits with Model 2 are given in Figure 5(b) The 



systematic effect of our choice of PnthComp will be dis 
cussed later by comparing with results using other values 
of TnthComp in the range of 2.0-2.7. 
We estimated Ah and E^ with the two Suzaku obser- 
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Fig. 6. — The ratios of the spectra on the two ends of each branch. 
Spectra with high total PCA intensity are divided by those with 
lower total PCA intensity in order to illustrate that on each branch 
the spectrum evolves differently. 

vations. The energy bands and the data bi nning method 
used were the same as in ILin et ail (|2010t ) . We applied 
Model 1 (Model 2 gives very similar results) and fitted 
both Suzaku observations simultaneously with Ah and 
Sga tied to be the same for both observations. Only a 
very weak CPL is needed in both observations. The ab- 
sorption edge cannot be constrained by the data, proba- 
bly because of the gap between the XIS and PIN (10- 
11 keV) and the relatively poor data quality around 
this edge. We still included it in the fit but fixed 
-Ecdgc = 8.93 keV (see above). Finally, we obtained 
A H = (2.35±0.01)xl0 22 cm- 2 and £ ga = 6.64±0.06 
keV. This A H is con sistent with A H = ( 2.38± 0.12) x 10 22 
cm -2 obtained by iWroblewski et all ([2008D from the 
modeling of absorption edges found in Chandra high- 
resolution X-ray spectra. 

4. RESULTS 
4.1. Spectral Evolution 

Figure [3] shows that the source goes back and forth 
between the HB and NB several times in the first seven 
days and then up and down the FB many times in the 
last two days, but no clear secular changes are seen in 
Figure [2] (upper panels). As in LRH09, we calculate the 
ratios of spectra in different branches, which are shown 
in Figure [Bl The blue dotted line is for the HB (the upper 
vertex divided by the top of the HB) and shows that the 
spectrum pivots around 9 keV, with the intensity below 
increasing and above decreasing as the source descends 
the HB. The green dashed line corresponds to the NB. 
Although the intensity increases over the whole energy 
range shown as the source evolves up the NB, the ef- 
fect becomes stronger above ~7 keV. The increase of the 
intensity on the FB is also mostly in the high energies 
10-30 keV (red dot-dashed line). The ratios for the NB 
and FB drop sharply around 30 keV, above which in fact 
there is no significant source emission. These ratios are 
similar to XTE J 170 1-462 in the Sco-like stage (LRH09). 

The parameters derived from the spectral fits are 
shown in Figures [?H9] and are listed in Tables [IHU For 
all luminosity (L) and radius (R) quantities, we assume 
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Fig. 7. — The fit results of Sz-resolved spectra. The luminosities of the thermal components are plotted against their characteristic 
temperatures. The dotted, dot-dashed, and dashed lines correspond to R = 8.0, 6.3, and 3.0 km, respectively, assuming Lx = 47r R? ctsbT 4 , 
and the red solid line is a constant M line (see text). Note that for Model 1 we plot the thermal disk emission unscattered and observed 
directly, while for Model 2 we plot the thermal disk prior to scattering. 



a source distance of 12.6 kpc (Section 14. 3|) . and the 
luminosities are based on bolometric fluxes of relevant 
spectral components (for the CPL component, we inte- 
grate over energies above 1 keV). Isotropic emission is 
assumed for all components (but see discussion in Sec- 
tion I5.2j) . except the MCD one, whose inclination i is 
assumed to be 60°. Rbb is the apparent radius of the 
BB component and is related to the normalization as 
N B b = (RBB±m/D wkpc ) 2 , where £>i k P c is the distance 
to the source in units of 10 kpc. Rmcd in Model 1 is 
the apparent inner disk radius and is calculated from the 
MCD normalization N M cu = (i?MCD,km/£>iok P c) 2 cos i, 
while i?MCDPS in Model 2 is from the MCDPS normal- 
ization A^MCDPS = (-RMCDPS.km/^lOkpc) 2 COS«. 

Figure [7] shows the luminosities of the thermal compo- 
nents MCD/MCDPS and BB versus their characteristic 
temperatures, while Figures [8] and [9] show the spectral 
fit results as a function of the rank number S%. We in- 
cluded several constant radius lines in Figure [3 assum- 
ing Lx = 47t<7sb^ 2 7 i4 - The dotted and dot-dashed lines 
have R = 8 km (iV BB = 40) and R = 6.3 km (iV BB = 25), 
corresponding to the apparent sizes of the NS from the 
sum of the boundary layer and the burst emission and 
from only the burst emission, respectively (Section I4.3[) . 
The dashed lines have R — 3 km, shown for reference. 
The red solid lines describe the relation between I/mcd 
and /cTmcd with varying i?MCD at a constant M, i.e., 

Lmcd oc T^ctj (LRH09; referring to the MCDPS in 
Model 2). 

We focus on the results of Model 1 for the HB first. 
Figures [THE] (left panels) show that Lmcd decreases and 
Lcpl increases as the source climbs up the HB, while 
L BB changes much less. The anti-correlation between 
Lmcb and Lcpl is clear, as their sum (filled blue squares 
in the left panel d of Figure [8]) changes by <10% on the 
HB, much less than either individual component. 



Model 2 assumes that some corona surrounds the ther- 
mal disk and produces the Comptonization emission. 
From this model, we also see the strong anti-correlation 
between the (unscattered) MCD component and the 
Comptonization emission (Figure [8]). However, unlike 
Model 1, we see that the apparent inner disk radius from 
Model 2 (inferred from the MCDPS) increases as the 
source climbs up the HB. Moreover, the thermal disk 
based on the MCDPS follows a constant M line (the red 
solid line in Figure [7|) • 

The Comptonized component becomes weak on the NB 
and is almost negligible on the FB (Figures [5] and |8} , 
and the results of other components become very similar 
between Models 1 and 2 in these two branches. In Fig- 
ures [7HHI both models show that in the NB the largest 
changes occur in the BB component. Its normalization 
varies by a factor of about 2.5, with the temperature hov- 
ering around 2.6 keV (Tables[l]-[2]). This suggests that the 
change of intensity above ^7 keV (green line in Figure [5]) 
is mostly due to the BB component. However, since the 
BB is just a small component (<30% in terms of lumi- 
nosity), the total luminosity only changes by ~20% on 
the NB (Figure 1 and Tables HH2J) ■ 

The variations of the MCD/MCDPS component on the 
FB approximately follow a constant M line (the red solid 
line in Figure[7]) from both models. The disk temperature 
increases from ~1.7 keV at the lower vertex to nearly 3 
keV at the top of the FB. Compared with the dotted line 
for the size of NS, the disk is consistent with being trun- 
cated outside the NS. The results of the BB component 
show large uncertainties, especially in the upper part of 
the FB, making it difficult to infer the evolution trend of 
this component in this branch. This was also seen in the 
spectral fits of XTE J 170 1-462 in LRH09. 

We have seen that the MCDPS from Model 2 evolves 
approximately along a constant M line over the whole Z 
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Fig. 8. — The fit results of Sz-resolved spectra as a function of the rank number Sz- The error bar of Sz corresponds to the range of 
Sz for each spectrum. Ltot is the total luminosity, whose error bar is not plotted but is given in Tables [TK2l See text for the meanings of 
other quantities. The units of the R, kT, L quantities are km, eV, and 10 erg 



track. The fractional rms of M using the relation of M oc 
^mcdps-Rmcdps is about 6%. Assuming Lmcdps = 
k^MCDPS> we estimate the slope by minimizing the \ 2 
value defined as J2 (Vi ~ P x i ~ a ) 2 / + P 2 (4i + 4 ) ) ( e -g-, 
iTremaine et al.l [2002T ) . where y = log(LMCDPs) with la 
error e y , x = log(T M CDPs) with la error e x , a = log(fc), 
and the subscript i denotes the data points used, eo is 
introduced to account for the possible systematic error r 
of /cTmcdps (eo = log(l + r)) and is set to such a value 
that the reduced \ 2 is one - We obtain j3 = 1.36 ± 0.07 



respectively. 

from data of the whole Z track, with r — 2.4% used. 
Consistent values of j3 can also be obtained using other 
choices of r nt hComp in the range of 2.0-2.7 in Model 2. 

Figure [9] shows the variations of the Fe K line and 
absorption edge using Model 1, while Model 2 gives very 
similar results. The line luminosity Lii ne first slightly 
decreases from the top of the HB to the upper vertex 
and then increases again until the end of the FB. It is 
highest on the FB. The line width au nc is ~0.4 keV. 
The equivalent width EWu nc peaks near the lower vertex, 
with a value of ~ 120 eV, and is ^60 eV around the upper 
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Fig. 9. — The fit results of Gaussian Fe line as a function of the 
rank number Sz using Model 1. 
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Fig. 10. — The fractional rms versus the luminosity fraction of the 
Comptonized component, i.e., nthComp, obtained with Model 2. 

vertex. The absorption depth at the edge energy T d g0 
is around 0.05 and peaks at the lower part of the FB. 
We note that as RXTE data have only a modest energy 
resolution (~1.2 keV FWHM at 6.5 keV) the results of 
the Fe K line and absorption edge should be regarded 
with caution. 
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Fig. 11. — The dependence of the frequencies of the upper kHz 
QPO and the HBO on S z . 
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Fig. 12. — The frequency of the upper kHz QPO versus the inner 
disk radius (prior to Compton scattering) from Model 2. The inner 
disk radius is derived from the MCDPS normalization, without 
corrections for factors such as the hardening effect. The solid, 
dotted, and dashed lines show the relation of frequency oc ^MCDPS' 
with /3 = —1.5 (corresponding to the Keplerian frequency at the 
inner disk radius), —1.72 (from the fit to the HB data), and —2.01 
(from the fit to the HB and NB data), respectively. 



4.2. Comparison with Fast Variability 

Figure HOI shows the rms integrated from the PDS (0.1- 
10 Hz) versus the Comptonization fraction inferred from 
Model 2. The same type of plot, substituting results 
from Model 1, appears very similar. Both models infer 
low Comptonization fraction (<10%) in the NB and FB, 
with significantly higher values for the HB. The rms vari- 
ations appear to be roughly correlated with the Comp- 
tonization fraction, and only the HB has rms >5%. In 
terms of these properties, the NB and FB are similar to 
the soft /ther mal state of bla ck hole X-ray binaries and 
atoll sources ([Lin et al. 2007!), while the HB is closer to 
their transitional/steep power- law state. We note that 



the rms shows an enhancement near the lower vertex, 
which is due to the normal branch oscillations. We also 
note that the FB in fact shows large variability in light 
curves. The low rms obtained above to some extent is 
because we only integrate the rms above 0.1 Hz. Using a 
lower integration limit could result in a larger rms, as the 
PDS in the FB is st eep, with a power-law slope around 2 
(jHoman et al.ll2002l ). I n the above, we j ust use the same 
frequency range as in ILin et all (12001 and LRH09 for 
fair comparison with their results. 

Table |3] lists the HBOs (fundamental) and upper kHz 
QPOs that we detected in our data set. The dependence 
of their frequencies on Sz is plotted in Figure [IT] There 
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Fig. 13. — The frequency of the HBO versus the apparent inner 
disk radius (prior to Compton scattering) from Model 2. The solid, 

dotted, and dashed lines show the relation of frequency oc R^Qj^pg, 
with /3 = —3 (corresponding to Lense-Thirring precession), —3.70 
(from the fit to the HB data), and —4.35 (from the fit to the HB 
and NB data), respectively. 

are also other QPOs such as HBO second harmonics and 
lower kHz QPOs, which are not explore d here due to 
their smaller ranges spanned in the Z track (jHoman et al.l 
12002ft . The HBOs span the whole HB and NB ranging 
from about 20 to 60 Hz, while the upper kHz QPOs span 
the whole HB and a small upper part of the NB rang- 
ing from abo ut 620 to 1000 Hz, ge nerally consistent with 
the results in lHoman et al.l (|2002T ) . However, they found 
some upper kHz QPOs from around 1000 up to 1090 Hz 
on the middle NB, which we could hardly detect from 
our data. This is probably because they used more ob- 
servations in their data analysis plus the fact that these 
QPOs are weak to begin with (rms < 2%). The HBOs 
become broader on the lower NB, and their FWHM val- 
ues have lar ge uncertainties. Th us we fixed the FWHM 
to be 15 Hz (jHoman et al.| [2002) in some fits. 

Figure IT21 shows the frequency of the upper kHz QPO 
versus -Rmcdps from Model 2. Assuming a relation of fre- 
quency = fc^MCDPS ano - using the slope estimate method 
in SectionEm we obtain f3 = -1.72t"' 3 4 l (the dotted line 
in Figure [12] a systematic error of 2.3% is added for 



-0.35 
-0.51 



-Rmcdps) using only data on the HB and f3 — —2.01 
(the dashed line; a systematic error of 2.1% is added for 
-Rmcdps) if data on the NB are also used. These values 
are consistent with —1.5 (the solid line, forced to pass 
the data point with Sz=0.96-1.00), expected if the fre- 
quency of the upper kHz QPO is the Keplerian frequency 
at -Rmcdps ■ The systematic error of f) due to our choice 
of r nt hComp in Model 2 is estimated to be around 0.6 
based on comparison with results using other values of 
rnthComp in the range of 2.0-2.7. In Model 1, -Rmcd on 
the HB increases with the frequency of the upper kHz 
QPO. Therefore, it is important to note that the rela- 
tionship between the upper kHz QPO frequency and the 
inner disk radius changes sign between the use of Mod- 
els 1 and 2. The result shown in Figure [T"2l depends on a 
model that counts the number of Comptonized photons, 



so that the QPO frequency becomes correlated to the 
apparent inner radius of the disk prior to the effects of 
Compton scattering. 

Figure [13] shows the frequency of the HBO versus 
-Rmcdps from Model 2. We see that the frequency of 
the HBO increases as -Rmcdps decreases down the HB, 
but both become relatively constant on the NB. We note 
that in fact the frequency of the HBO decreases slightly 
on the lower NB (Figure ITT1 iHoman et all 12002ft . Also 
assuming a relation of frequency oc R 



MCDPS' 



we esti- 



mate P = -3.701^ (the dotted line in Figure [£3] a 
systematic error of 1.6% added for -Rmcdps) using only 
data on the HB and (3 — — 4.35+ ; 94 (the dashed line; a 
systematic error of 2.8% added for -Rmcdps) if data on 
the NB are also used. These values are roughly consis- 
tent with —3 (the solid line in Figure [IB"]) , a value ex- 
pected in a Lense-Thirring interpretati on of the HBO 
([Stella fc VieTrl[T99l Ivan der Klisl 12006ft . Based on the 
results using other choices of r nt hCom P , we estimate a 
systematic error of 1.2 for the above measurements of /3. 

4.3. The Source Luminosity Based on Type I X-ray 

Bursts 

Type I X-ray bursts in NS LMXBs can provide criti- 
cal reference quantities such as the Eddington flux and 
the apparent NS surface area when they show photo- 
spheric radius expansion. These quantities can be com- 
pared with the spectral modeling results to infer the 
source luminosity and the visible size of the boundary 
layers, which are key elements in c onstructing the physi- 
cal picture of the accretion process (|Lin et al.l l2007. 2010, 
LRH09). GX 17 +2 shows both short (~10 s) and long: 
(>10 s) bursts (jKuulkers et al.l 120021 : iGallowav et~al1 
2008). With the persistent emission comparable to the 
peak net burst emission, these bursts are probably ob- 
served at accretio n rates around the Eddington limit and 
are n ot expected (|Gallowav et al J 120081 : iRemillard et al.l 
2006). It is still under debate on how to analyze bursts 
at such high accretion rates. In the "standard" burst 
analysis, burst spectra are obtained by subtracting the 
persistent emission from the total emission and then fit- 
ted with a simple BB model. At high accretion rates, 
the persistent emission from the boundary layer may be 
strong. If it varies during bursts, it is instead proba- 
bly more appropriate to fit the total emission with mul- 
tiple spectral components, typic ally including a BB to 
d escribe the burs t emiss ion (e.g.. iSztajno et aLlll986ft . 

Kuulkcr s et al.l (|2002l ) applied the "standard" analy- 
sis to ten bursts of GX 17+2 from RXTE data and ob- 
tained successful fits for most of them. Here we combine 
our spectral modeling results of the persistent emission 
with their burst analysis to infer the Eddington flux and 
apparent NS area. In our data set, there are no short 
b ursts but five long one s, corresponding to bursts 6-10 
in lKuulkers et all (|2002l ). Burst 10 is near the lower ver- 
tex (Sz = 2.18) but shows no photospheric radius expan- 
sion. Bursts 6-9 are near the upper vertex (Sz = 0.86, 
1.04, and 0.73, respectively) and show photospheric ra- 
dius expansion, and we focus on them. In these bursts, 
for several tens of seconds right after the touch-down 
(the moment when the expanded matter falls back onto 
the NS surface), the burst temperatures peak and hover 
around ~2.65 keV, and the apparent net burst area is 
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roughly constant (Nbr— 25), corresponding to a peak net 
burst flux around 1.24xl0 -8 erg cm -2 s _1 (from burst 
6). In comparison, the boundary layer during the per- 
sistent emission has a similar temperature and Abb— 15 
from our modeling, about a size of 60% of the net burst 
area. 

With such a strong boundary layer, the Eddington flux 
and the apparent NS area would be very different de- 
pending on whether we infer them only from the net 
burst emission or have the boundary layer emission also 
included. We note this is not a serious issue in atoll 
sources, as their boundary layers are small (<10% of the 
NS surface) and have temperatur es typically lower than 
the bursts (jLin et al.ll2007l I20T0I LRH09). Considering 
that both the boundary layer and the burst emission are 
from the NS surface, we add them together and obtain 
the Eddington flux of about 2.0xl0 -8 erg cm -2 s _1 and 
the total apparent NS surface area of Abb— 40 for GX 
17+2, but we also discuss the results with the Edding- 
ton flux and the total NS apparent area calculated from 
the burst emission only. 

The luminosities in Tables [TH5] are in units of the Ed- 
dington flux obtained above, and they would be 60% 
larger if we do not include the boundary layer emission in 
the calculation of the Eddington flux above. In the upper 
vertex, GX 17+2 has luminosity of a factor of 1.5 of the 
Eddington limit, higher than XTE J1701-462 in the Sco- 
like stage (below 0.9 Eddington luminosity (LRH09)). If 

M accounts for different source types (Section [5. ip . one 
might expect their luminosities to be similar. The above 
difference might be due to several uncertain factors, such 
as the inclination and the thickness of the accretion flow 
not taken into account in the Eddington flux correction 
above. However, we cannot rule out that their difference 
in luminosity is real , which could be due to factors such 
as different compositions of the accreted matterials and 
different NS masses. 

We estimate the source distance based on the Ed- 
dington flux obtained above. Using the empirical value 
of the Eddingt o n lum inosity of 3.79xl0 38 erg s" 1 from 
iKuulkers et all pOOl . the source distance is 12.6 kpc, 
which is assumed in this paper. If the theoretical 
expression of the E ddington limit (see Equation 8 in 
IGallowav et al.|[2008f) is used and a NS with a mass of 1.4 
M & and a radius of 10 km is assumed, we obtain a source 
distance of 10.5 kpc for the H-poor case (the H-fraction 
X = 0) and 8.1 kpc for the H-rich case (X = 0.7). The 
above distances should increase by 26% if we only use 
the net burst emission to estimate the Eddington flux. 

5. DISCUSSION 

5.1. The Role of M in the Spectral Evolution 

The spectral evolution o f GX 17+2 along the Z track 
based on Model 1 (Figure |7(a)[ ) is very similar to that 
of XTE J1701-462 in the ScT>Tike stage as obtained by 
LRH09 using a similar model. This supports the viability 
of physical interpretations for Z branch evolution offered 
in that paper (see Section Hx2"|) . An associated question 
is whether the M into the disk varies along the Z track. 
LRH09 inferred a constant M along the Sco-like Z tracks 
of XTE J1701-462 based on the MCD component for the 
NB and FB. Although in the HB the MCD component 
did not follow a constant M, the variation of the MCD 



luminosity was largely offset by that of the Comptonized 
component, leading LRH09 to surmise that on the HB 
M is also constant, with the thermal disk emission sim- 
ply converted to the Comptonized emission as the source 
ascends the HB. Similar behavior observed in GX 17+2 
using Model 1 implies that M is probably also constant 
in this system. Model 2 took a step further by providing 
an empirically self-consistent way to model the Comp- 
tonized component, assuming that a hot corona in the 
line of sight to the disk Comptonizes part of the thermal 
disk emission. We have seen an intriguing result of this 
model, i.e., £mcdps, the thermal disk luminosity prior to 
Compton scattering, is consistent with lying along a con- 
stant M line over the whole Z track (Figure [7]). We note 
that in the above we use the MCDPS to infer the M into 
the disk. This assumes no significant mass loss around 
the area where most of the thermal disk emission is pro- 
duced. At the very inner part of the disk, some of this M 
can be taken away by mass outflow (Section I5.2[) . with 
the rest going onto the NS surface to form the boundary 
layer. Based on the littl e variation of t h e tota l flux on 
the HB and NB seen in IDi Salvo et~aT1 (120001) and the 
comparison with black hole X-ray binaries. iHoman et al.l 
( 2002) also suggested a constant M along the Z track. 

The conclusion of a constant M into the disk along 
the Z track is further supported by the global evolution 
of GX 17+2. Its secular changes were observed to be 
very small, either from our data set spanning nine days 
or f rom data spanning years in IWiinands et aT] (1997) 
and IHoman eTaTl ([200l . This is in contrast with strong 
secular changes seen in XTE J 170 1-462 in 2006-2007, 
which were a scribed to the variation in the M into the 
disk (LRH09. lHonran et al .1 [2007L 1201(1 . Then a reason- 
able explanation for the lack of strong secular changes in 
GX 17+2 is that its M has been fairly constant. This in- 
terpretation implies a need of hardly varying M over its 
entire Z track, in accordance with our spectral fit results. 

5.2. The Interpretation of Z Branches 

If the M into the disk does not change along the Z 
track as argued above, the remaining main question is 
what causes the evolution along the Z track. Based on 
XTE J1701-462, LRH09 suggested that the three Sco-like 
branches are due to three physical processes operating at 
a constant M into the disk. The source climbs up the HB 
due to the increasing portion of the thermal disk emission 
Compton scattered. The excursion from the lower to the 
upper vertices on the NB corresponds to the tran sition 
from a slim-disk accretion (lAbramowicz et aLlll988l) to a 
standard-thin-disk accretion (jShakura &: Sunvaevll 19731) . 
At the lower vertex, the disk is truncated at a radius out- 
side the ISCO as the local Eddington limit is reached at 
the inner disk radius (see discussion below). The source 
climbs up the FB due to the temporarily fast decrease of 
the inner disk radius toward the ISCO. Since the source 
evolved slower and stayed longer in the two vertices than 
in other parts of Z tracks (at least on the NB and FB), 
the above three processes could correspond to three kinds 
of instability (LRH09). The slim disk solution was asso- 
ciated with the upper vertex because it can provide more 
mass inflow onto the NS surface than the standard thin 
disk and can explain the increase of the BB apparent area 
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Top of the HB 




Fig. 14. — Schematic sketch of the accretion flow at different 
positions of the Z track. The system is assumed to be axisymmetric 
relative to the rotational axis, and the sketch shows one half of the 
cross-section through the rotational axis. Arrows represent the 
mass flow onto the NS surface and out of the disk (at the inner 
disk), with the strength roughly indicated by the thickness of the 
arrows. The blue arc represents the boundary layer (BL). 

from the lower to the upper vertices. The disk emission 
was not observed to change significantly during such a 
transition (in the Sco-like tracks), probably because the 
slim disk emission differs from the standard disk signifi- 
cantly only at accretion rates much higher than the Ed- 
dington limit ()Mineshige et alj|2000h . 

As a similar spectral evolution is observed in GX 17+2, 
the above picture can be applied to this system and is 
sketched in Figure [14l The disk is drawn to be thicker for 
the HB to represent a slim disk, and is thinner for the FB 
to represent a standard thin disk. On the HB, the disk 
is surrounded by a Comptonizing corona with a temper- 
ature of about 6 keV, based on our fits with Model 2 in 
this work. Nearly one half of the thermal disk is scat- 
tered at the top of the HB, but this quantity has a strong 
dependence on r nt hComp (Section [3]). The optical depth 
of the corona is about a few from Model 2, based on its 
relati on with r nt hCom P and fcT ei „thCom P (|Zdziarski et al.l 
1996). The Comptonized component is strong only at 
the top of the HB, making it hard to measure the vari- 
ation of the optical depth along the Z track. Thus we 
cannot easily infer whether the decrease of the Comp- 
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tonization emission as the source descends the HB is due 
to the decrease of the covering fraction of the corona or 
due to the decrease of its optical depth. 

We add some small arrows in Figure to indicate the 
possible strong hot mass outflow at the very inner part of 
the accretion disk in this system. This is expected to be a 
common phenomenon for a disk accretion system at high 
accretion rates, in which the mass can be input into the 
disk from outside and ejected as a disk wind at the very 
inner part of the disk where the emission perpend icular 
to the disk plan reaches the local Eddington limi t (IKata 
1980 | IWatarai et al.l [200CH IMineshige et al.| [2000t iFukuel 
2004t lOhsuga fc Mineshigell2007l) . The mass outflow can 
also explain the low BB luminosity, relative to the disk 
luminosity (discussed below). 

Figure [TJ] does not include the jet component that 
is possibly responsi ble for the rad i o em ission detected 
from this system. IMigliari et al.l (]2007f ) observed the 
radio emission to decrease significantly from the HB 
to the FB along the Z track and found it to correlate 
with a hard spectral component (the "hard tail"), which 
they fitted with a PL. In this paper we used a differ- 
ent dat aset and a diff e rent spectral model from those 
used by [Migliari et al. (2007). However, we found that 
consistent results can also be obtained if we fitted our 
data with the i r mod el (the same as the one used by 
iDi Salvo et "all (120001 ) V The CPL/nthComp component 
in our model was also observed to decrease significantly 
from the HB to the FB, i n a similar way as the hard 
tail in Mi gliari et al.l (120071 ) . thus probably also correlat- 
ing with the radio e mission. It played a s imilar role as 
the PL component in IMigliari et al.l (|2007l ) in that both 
of them are the dominating spectral component at en- 
ergies above 30 keV when they are the strongest, at the 
top of the HB. We have shown that the CPL/nthComp 
model described our data better than the PL model (Fig- 
ure 01) • We note that the HB spectra in IMigliari et al.l 
(2007) seemed to extend above 60 keV, which was not 
seen in our HB spectra. The 60-100 keV HEXTE count 
rates are about three a and about 2% of the background 
level in their HB observation 70023-01-01-00. To check 
whether the apparent lack of the hard tail above 60 keV 
in our data is due to relatively small exposures of our Sz- 
resolved spectra (each with a few ks at the top of the HB, 
compa red to ~10 ks for the spectrum in Migliari et al. 
(|2007h ). we combined the spectra of the first five Sz se- 
lections (representing the upper half of the HB, where 
the hard tail was reported to be strongest), resulting in 
a spectrum with an exposure time of ~20 ks. We found 
that the source was detected at less than one a in the 
60-100 keV band (combining data from both HEXTE 
clusters), suggesting that this (additional) hard emission 
may be transient, even for the same location in the Z 
track. 

The hard tail observed by BeppoSAX from GX 17+2 in 
the HB has been explained in two different models, i.e., 
hybrid Comptonization and bulk motion Comptoniza - 
tion, bv lFarinelli et al.l (|2005| ) and iFarinelli et al.l (|2007f ). 
respectively. In the hy brid Compton ization model (us- 
ing the eqpair model (|Coppil 11999ft in XSPEC), two 
electron populations (thermal and non-thermal) are re- 
quired, with the non-thermal one accounting for most of 
the hard tail above ~30 keV. In the bul k motion Comp- 
tonization model (using the bmc model (|Titarchuk et al. 
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11997ft in XSPEC), the hard tail is caused by Comp- 
ton upscattering of soft photons by energetic electrons 
in a converging flow onto the NS. The hard component 
(CPL/nthComp) in our HB spectra has an exponential 
cutoff at ~10 keV and has been explained by us as due to 
Comptonization by a low-temperature (~6 keV) thermal 
plasma covering the accretion disk (i.e., the non-thermal 
electron population is not needed). The bulk motion 
Comptonization model cannot explain this component 
easily because the cutoff energy of its emergent spec- 
trum is much higher (a few hundred keV. lTitarchuk et al.l 
IT997I) . 

In Figure Q3] the boundary layer, which we model with 
a BB, is pictured as an equatorial belt on the NS sur- 
face. On the HB, Nbb of the boundary layer is around 
15, while the entire NS surface has a corresponding value 
of 40 if we sum the boundary layer area and the burst 
emission area together (Section 14. 3p . indicating that the 
boundary layer has an apparent fraction of 3/8 of the 
entire NS surface. The real fraction is different, because 
of the special geometry of the bound ary layer picture d 
above, and can be estimated following ILin et aLl {2007). 
We infer the maximal latitude of the boundary layer to 
be 26° and 35°, corresponding to 43% and 57% of the 
entire NS surface, for an inclination of 60° and 30°, re- 
spectively. In the lower vertex, the boundary layer has 
^Vbb~6, and we infer the maximal latitude of the bound- 
ary layer to be 11° and 18°, corresponding to 19% and 
30% of the entire NS surface, for an inclination of 60° 
and 30°, respectively. The solutions for an inclination of 
60° are used to create Figure Q2J except for the panel for 
the top of the FB, which is arbitrary due to large uncer- 
tainty. The above estimates of the size of the boundary 
layer are only approximate, neglecting factors such as the 
thickness of accretion flow, and should be larger if we es- 
timate the NS surface area only from the burst emission 
area. The above corrections imply that the real bound- 
ary layer luminosities should be higher than those shown 
in Figures [THE] by up to 50%. However, even with such 
corrections, the luminosity of the boundary layer seems 
to be still less than the disk luminosity, which is not ex- 
pecte d from the simple energy argument (Mitsuda et al. 
1984). One explanation for this is the mass outflow de- 
scribed above. 

Although large uncertainties are seen on the FB (Sec- 
tion 14. 1|) , the boundary layer on the HB and NB has 
temperatures around 2.6 keV, about the peak value 
seen in the radius ex pansion bursts from this system 
(|Kuulkers et al.ll2002ft . This indicates that the bound- 
ary layer emission may have reached the local Eddington 
limit. In such a situation, the increase in the mass ac- 
cretion rate onto the NS surface will just lead to increase 
in the boundary layer emission area (not its tempera- 
ture), as seen on the NB. This was also observed in XTE 
J1701-462 (LRH09). 

5.3. The Origins of kHz QPOs and HBOs 

We have found that the dependence of the frequency 

of the upper kHz QPO on the apparent inner disk radius 

3/2 

-Rmcdps is roughly consistent with frequency oc -R MCDPS 
from Model 2, making it natural to identify it as the Ke- 
plerian frequency at the inner disk radius. This has been 
a critical assumption in many explanations of the kHz 



QPOs (|Lin et al.ll201~H Ivan der Klidl200ll . We note that 
-Rmcdps has been corrected for the photons which we as- 
sumed to be from the thermal disk emission but scattered 
by the corona. We also caution that the absolute values 
of the radius measurements have scale uncertainties that 
depend on the source distance , the inclination, th e hard- 
ening effect, and the NS mass (iZhang et al.lU997l). Th ere 
is another important factor (rj in IZhang et al" [ Tl99l to 
account for the difference between the real inner disk ra- 
dius and the radius where the disk temperature peaks. 
All these factors have too large uncertainties to allow for 
a meaningful estimate of the real inner disk radius di- 
rectly from .Rmcdps- We calculated the variation of the 
upper kHz QPO with -Rmcdps to infer its dependence 
on the real inner disk radius, assuming that the above 
conversion factors bewteen -Rmcdps and the real inner 
disk radius are fairly constant, at least along the Z track 
where upper kHz QPOs are observed. 

The frequency of the HBO is too low to be the Keple- 
rian frequency at the inner disk radius. However, we also 
observe its close dependence on -Rmcdps covering the en- 
tire HB and NB. Thus the HBO is probably intimately 
related to the dynamics at the inner disk radius in some 
way. In the Lense-Thirring interpretation for the HBO, 
an app roximate relation of frequency oc flunn pg is ex- 
pected (|Stella fc Vietrilll998tlvaii der Klidl2009 ). and our 
results are roughly consistent with this. Although the 
frequency of the HBO varies much less on the NB than 
on the HB, it still shows a small increase first and then 
a small decreases as the s ource descends the NB (Fig- 
ure [TTJ iHoman et al.l l2002h . The error bars of -Rmcdps 
are not small enough for us to conclude whether it cor- 
relates with this evolution pattern. However, we cannot 
rule out that such an evolution of the HBO on the NB 
is due to the change of the disk structure from a slim 
disk to a standard thin disk on the NB (Section 5.3). 
The decrease in the quality of the HBO on the lower NB 
(Ho man et al.ll2002|) might be related to this. 

6. CONCLUSIONS 

We have analyzed a Z track of GX 17+2 from RXTE 
observations between 1999 October 3-12. Our analysis 
has confirmed similar spectral properties of GX 17+2 to 
those of XTE J1701-462 in the Sco-like stage. In this 
study, we took a further step by modeling the Comp- 
tonized component in an empirically self-consistent way 
and found that the evolution of the thermal disk is con- 
sistent with being at a constant M into it over the whole 
track, supporting our conclusion of a constant M over 
the whole Sco-like Z track. Similar to XTE J1701-462, 
we found that the branches of GX 17+2 can be explained 
as being due to three processes operating at a constant 
M into the disk. The HB is due to increase of the Comp- 
tonization with respect to the upper vertex. The FB is 
due to shrinking the inner disk radius on fast timescales 
from the lower vertex, on the NB we saw variation in the 
area of the boundary layer, which can be explained as be- 
ing due to transition from accretion through a slim disk 
in the upper vertex to a standard thin disk in the lower 
vertex. The boundary layer has a size of ^20-30% of the 
entire NS surface in the lower vertex and ~40-60% in 
the upper vertex and is probably at the local Eddington 
limit. 
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The dependence of the frequency of the upper kHz 
QPO on the apparent inner disk radius is found to be 

roughly consistent with frequency oc R 



-3/2 
MCDPS 



on the HB, 



supporting identification of the frequency of the upper 
kHz QPO as the Keplerian frequency at the inner disk 
radius. We measure frequency oc #mcdps' 7 f° r tnc HBO 



over the entire HB and NB, indicating an intimate re- 
lation between the HBO and the dynamics at the inner 
disk. 
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TABLE 1 
Spectral fit results using Model 1 
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Note. — The first column is the Sz range of each spectrum. The second column is the exposure per PCU. The last column is the total luminosity 
in units of the Eddington luminosity (3.79xl0 3S erg s — 1 ; Section |4.3[ l . Sec text for the meanings of other columns. 
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TABLE 2 
Spectral fit results using Model 2 
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0.00+° 


05 


0.40±0 


12 


2.2±0 


4 


84±19 


6.3±1 


2 


1 


33 


118) 


1 


97 ± 0.10 


3.80- 


-4 


00 


2.39±0.12 


49+1 




2.91±0.27 


6 n+ 6 - 3 

O.U_ 3 2 


0.04+° 


04 


0.40±0 


13 


2.5±0 


5 


88±19 


6.6±1 


3 


1 


24 


128) 


2 


06 ± 0.16 


4.00- 


-4 


20 


2 54+0.O8 


42+_l 


43+3 


3.16±0.35 


Q r + 5.8 
J -°-0.5 


0.01+° 


07 


0.41±0 


13 


2.4±0 


4 


83±20 


5.3±1 


3 


1 


38 


121) 


2 


14 ± 0.15 


4.20- 


-4 


40 


2 62+ " 11 


3S +_l 




3 14+0-00 


3 fi+ 8 - 2 


0.02+° 


06 


0.44±0 


14 


2.6±0 


5 


83±20 


5.4±1 


3 


1 


12 


125) 


2 


23 ± 0.21 


4.40- 


-4 


80 


2 89+ - 05 

^•o»_ 0.31 


29±l 


291;° 


3.49±0.71 


1 1+9-3 
i ' J -_ 0.2 


0.00+° 


05 


0.44±0 


13 


2.5±0 


5 


77±21 


4.2±1 


2 


1 


02 


127) 


2 


33 ± 0.24 



Note. — Sec text and Tablc[^for the meaning of each column. 



17 



HBO fundamental 
Sz Frcqucncy(Hz) FWHM(Hz) Fn^i" 



0.07- 





.24 


22.0±0.3 


4.3±1.3 


0.24- 


-0 


.32 


24.2±0.5 


4 8+?I 

— 1.6 


0.32- 


() 


.40 


26.9±0.3 


5.7±1.2 


0.40- 


(j 


.48 




o r + 13 
J -°-0.9 


0.48- 





.56 


29.8±0.0 


4.4i 1 .6 


0.56- 


-0 


64 


34.2±0.7 


4 3+5-9 

*-°-2.3 


0.64- 





.72 


o6.u±0.4 


7.4il.6 


0.72- 


-0 


.80 


39.6±0.5 


8 7+?! 


0.80- 


() 


.88 


42.3±0.5 


9.4±2.1 


0.88- 


-0 


.96 


45.1±0.7 


9.2±3.0 


0.96- 


-1 


.00 


49.6^;° 9 




1.00- 


-1 


.04 


52.3±0.9 


"t" 


1.04- 


-1 


12 


53.3±1.7 


o 1+12.8 
'■'-1.9 


1.12- 


-1 


20 


57.2±0.9 


5.9±2.7 


1.20- 


-1 


.28 


59.6±0.7 


3.3±2.4 


1.28- 


-1 


.36 


60.0±0.7 


4.0±3.2 


1.36- 


-1 


.44 


60.9±1.8 


g+ 19 - 5 


1.44- 


-1 


.52 


56.4±5.7 


15.0 


1.52- 


-1 


.60 


53.3±4.0 


16 4+ 117 


1.60- 


-1 


.68 


55.3^ 9 x 


15.0 


1.68- 


-1 


.76 


54.51H 


13.7t^ 7 


1.76- 


-1 


.84 


53.8+ 7 ° 2 5 


15.0 


1.84- 


-1 


.92 


51-o+t; 


15.0 


1.92- 


2 


.00 


44 5+ 5 - 


15.0 



TABLE 3 

HBO FUNDAMENTAL AND UPPER KHZ QPOs 



Lower kHz QPO 

Froquoncy(Hz) FWHM(Hz) rins 



Upper kHz QPO 

Frequency(Hz) FWHM(Hz) rins 



0.032±0.005 

0.029±0.006 

0.035±0.004 

0.024±0.004 

0.028±0.005 

020+ - 013 
u.uzu_ 005 

0.030±0.003 

0.028±0.003 

o n27 +0 - 003 

U.U^i ( _o,005 

0.023±0.004 
014+ - 013 

U ' UA ^— 0.006 

0.017±0.005 

009+ ' 007 
u.uu»_ 004 

0.015±0.003 
0.014±0.003 
n m a+o.004 
0.018_ 007 

0.011±0.004 
0.017±0.004 
0.023±0.005 
0.018±0.005 
0.027±0.008 
0.022±0.006 
0.020±0.006 
0.032±0.007 



,+8.1 

552.6±2o"3 



518.2! 



619.5±22.9 
653.6±12.0 
690.3±16.1 



,+21.5 
-12.0 

58.1±36.0 



22.1] 



67.1±35.9 

48 3+38-8 

70.9±28.6 



0.015±0.004 
0.017±0.004 

0.026±0.007 
0.021±0.005 
0.024±0.004 



622.8±28.3 
626.0±21.1 

683.1± 20 " 4 6 
690.5±17.3 
693.0±32.4 
741.9±13.4 
773.6±23.8 
797.0±8.4 
827.3±12.0 
853.6±12.9 
899.3±26.2 
891. 3±\ 9 7 7 6 
962.7±45.5 
947.1±16.0 



232.3l« 7 ; 7 
149.5™ 

ii2.4±S;» 

156.9±50.1 

143 g+ 1112 

79 3+64 9 
' y '°-28.9 

205.5±85.0 
108.9±25.0 
123.1±35.0 
110.9±35.8 
141 4+ 104 - 4 

111-1^65.9 
34 fi+ 188 - 2 

174.5±110.5 
65.0t«; 7 



0.064±0.012 
0.055±0.010 
0.042±0.007 
0.049±0.007 
0.042±0.010 
0.039±0.008 
0.047±0.008 
0.040±0.004 
0.036±0.005 
0.032±0.005 
0.036±0.010 

n m s+o.016 
0.018_ 010 

0.028±0.013 

0.021±0.006 



